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Abstract 

We investigate the role of the effective nucleon-nucleon interaction in the description of giant 
dipole resonances in hot nuclei. For this purpose we calculate the response function of hot nuclear 
matter to a small isovector external perturbation using various effective Skyrme interactions. We 
find that for Skyrme forces with an effective mass close to unity an undamped zero sound mode 
occurs at zero temperature. This mode gives rise in finite nuclei (calculated via the Steinwedel- 
Jenssen model) to a resonance whose energy agrees with the observed value. We find that zero 
sound disappears at a temperature of a few MeV, leaving only a broad peak in the dipole strength. 
For Skyrme forces with a small value of the effective mass (0.4- 0.5), there is no zero sound at zero 
temperature but only a weak peak located too high in energy. The strength distribution in this 
case is nearly independent of temperature and shows small collective effects. The relevance of these 
results for the saturation of photon multiplicities observed in recent experiments is pointed out. 
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I. Introduction 



Giant dipole resonances built on nuclear excited states have been the subject of numerous studies 
since their first observation at the Berkeley 88" cyclotron in 1981 jl]]. By now a significant amount of 
information is available concerning their evolution with increasing excitation energy (which is mainly 
stored in rotational and thermal degrees of freedom) & ||, [|, |5| . The energy of the resonance is known 
to a good accuracy to be nearly independent of the excitation energy (E*) while the width increases 
up to approximatly E* ~130 MeV in the case of nuclei with A ~110. Beyond this value there are 
indications that the width saturates ||, ||, [| although some experiments suggest in contrast a steady 
increase ffl, H. In the recent experiments performed at the GANIL facility with a 37 MeV per nucleon 
Argon beam, the saturation of the width was observed at excitation energies greater than about 250 
MeV H together with a saturation of the photon multiplicity. At such high excitation energies, it is 
not clear, however, whether or not thermalization is actually reached. 

There are experimental indications that angular momentum contributes significantly to the width 
of the resonance up to the point where scission occurs @. Although some calculations of this effect 



1C, 11] show a rapid and continuous increase of the width with angular momentum, they do not 



seem to account neither for the observed variation at mean and high excitation energies nor for the 
saturation. 

At zero temperature the dominant contribution to the width is the spreading over two particle-two 
hole configurations [0|. It shows however a negligible variation with temperature [jl3L O]. In contrast, 



in Ref. 15], a rapid increase of the width at small temperatures and a saturation beyond T = 4 MeV 
was obtained by including explicitely those additional configurations (particle-particle and hole-hole 
states) which begin to occur in Random Phase Approximation (RPA) calculations as temperature 
rises. 

RPA calculations generally give resonance energies in agreement with each other (and with ex- 
perimental data). In contrast, rather different predictions can be found in the literature about the 
width dependence on temperature. Some RPA results for the giant dipole resonance for 40 Ca predict 
a sharp peak at zero temperature and an increase of the width with temperature |l(| whereas such is 



not the case for the self-consistent RPA calculations of Ref. [17]. In this last case, one already has at 
zero temperature a fragmented resonance displaying little evolution with temperature. In the nuclear 
matter calculations of Refs. [fill, O] it was found that for specific values of the particle-hole interaction 
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strength a disappearance of the collective (zero sound type) mode occurs at a temperature of a few 
MeV. This result was however obtained in the case of a schematic Skyrme-type interaction which 
contains no momentum dependent terms and therefore cannot be considered as a reliable effective 
force. 

The purpose of the present paper is to investigate what can be learned from nuclear matter calcu- 
lations of the isovector response function at finite temperature using more appropriate effective forces. 
In what follows we consider the standard effective Skyrme interactions which have been successful 
at describing accurately nuclear ground state properties while retaining contact with the standard 
fundamental description of nuclear matter based on the Brueckner reaction matrix ^(|. Calculations 
with such Skyrme forces are somewhat tedious. Nevertheless we will show that analytical formulae 
can be obtained which are still quite transparent. These formulae generalize those of Garcia- Reccio 



et al. [21] for the zero temperature response function. These formulae will allow us to discuss in 
rather general terms the evolution of the collective behaviour of nuclear matter with temperature. In 
particular we will see that zero sound is not as robust in nuclear matter as in usual Fermi liquids 
and disappears at temperatures of a few MeV for all the standard Skyrme forces investigated in the 
present work. 



II. The response function of hot nuclear matter 

We calculate the response function of nuclear matter at finite temperature to an infinitesimal external 
field of the form 

V ext = ET Z e-^ r e-^ + ^\ (2.1) 

where T3 is the third isospin Pauli matrix and r/ a vanishingly small positive number corresponding to 
an adiabatic switching of the field from the time t = —00. The temporal evolution of the one-body 
density matrix p is determined by the time-dependent Hartree- Fock equation in the presence of the 
perturbation term, i.e., 

ihd t p=[W + V ext ,p], (2.2) 

We consider effective Skyrme forces (see appendix A for details concerning the Skyrme force). In this 
case, the energy density can be written in terms of the one body density, kinetic energy density r and 
momentum density j . 
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A. Construction of the Response Function 

For a small enough external field it is legitimate to linearize the mean field evolution equation fl2.2j ) 
around its static solution. This procedure leads to the following approximate equation in momentum 
space for the difference 5p = p n — p p between the neutron and proton density matrices: 

ihd t (k\6p\k>) = (e(k) - e(k'))(k|Mk'} + 2(/(k') - /(k))(k|(W n - W P )\U) 
+4e(/(k') - /(k)) 5(k' - k - q) e -^+^)*. 

In this equation e(k) = h 2 k 2 /2m* (with m* being the effective mass; see appendix A) is the kinetic 
energy of a single particle state with momentum k in symmetric nuclear matter and /(k) is the 
corresponding occupation number 

/(k) = l/(l + e /3 W k )-^), (2.4) 

where (5 = 1/T is the inverse temperature. Note that since e(k) contains the kinetic energy only, the 
quantity p in the previous equation is not quite the chemical potential p but rather p — U, where 
U is the mean field defined in appendix A. From here on we will adopt the standard units in which 
h = c=l. 



The external field (2A) induces a difference between neutron and proton density distributions for 
which we consider the following time-dependent Ansatze suggested by the form of V eyi t 

(r\6p\r) = ae -il-* e -i(u+in)t 

{y\5t\y) = f3 e -ii-r e -i("+in)t (2.5) 

<r|«5J|r> = 7qe -iq-r e -i(o;+ir ) )t 
From the definitions of the densities given in appendix A the coefficients a, (3 and 7 must satisfy 

(a, A 7) = / (l,k.(k + q),l(2k + q).q)<k|M* = 0)|k + q> (2.6) 

The variation in the energy density can be written in terms of the quantities ( |2.5D , yielding 

W n (t) - W p (t) = 2V 5p(r, t) + 2V 1 X7.5p(r, t)V + 2V 1 5t(t, t) + 2W X (V.<Jj + tfj.V) (2.7) 

In this formula Vo and V\ are related to the parameters of the Skyrme interaction via the following 
expressions (see appendix A) 

^0 = ~% (xo + l)-k(x 3 + I) P% - ^(3ti(l + 2x0 + ta(l + 2x 2 )), 



2 V 2 J 12 V 2) ' u 16 v " ' (2. 

Vx = ^-(t 2 (l + 2x 2 ) -t 1 (l + 2x 1 )) 
16 



where po is the equilibrium density of nuclear matter. 

The retarded response function is determined by the corresponding polarizability, i.e., the ratio of 
the density change to the field strength 



II(w,q) = a/e. (2.9) 



In the case of the Skyrme effective force, this function is found by noting that Eq. ( |2.3| ) is solved by 
an Ansatz of the form 

6p(t) = 5p(t = Q) e -i("+iv)t } (2.10) 
provided the matrix elements at time zero satisfy 

(k\dp(t = 0)|k + q) = 4 jr [ { \ + q ; ' ~ ffl- . (V a + ^k.(k + q)a + e + Vrf - V l7 (2k + q).q) 

e(k + q) - e(k) + to + irj 

(2.11) 

By multiplying the previous equation respectively by 1, k.(k + q) and (2k + q).q and by integrating 
over k we obtain the following set of linear equations for a, (3 and 7 

a = (V + F1II2) a + T4/3II0 + 2m*Win 7 + ello 

P = (V U 2 + V1U4) a + V1U2P + 2 7 mWin 2 (2.12) 
2m* uj 

= 01 

1 q 2 (l- 2V 1 m*p ) 

The last equation for 7 can be checked to be a mere consequense of the equation of continuity for 
neutrons: 

OtPn + - V.tfj = l(h + t 2 ) ( Pn j - pj n ) (2.13) 

m 2 
with an analogous equation for protons. 

Solving the above linear system for a leads to the following expression of the retarded response 
function: 

IL(u, q) = _ n ° (u; ; q) 9 9 , g . (2.14) 

In this equation IIo is the unperturbed response function (often referred to as the Lindhard function 



[22]). In what follows the quantities II2 and H4 will be referred to as generalized Lindhard functions. 



They are defined as 



q) = 7^3 / d 3 fc J (k+ ^J^\ v (k-(k + q)) w , (2.15) 
(27r) d J a; + 177 - e(k) + e(k + q) 
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where the limit rj — ► + is implicit. Analytical expressions for the real and imaginary parts of T12n(u, q) 
are given in appendix B as well as the relation between our definitions and those of other authors plf , 
InEq. ( [HP we have used a modified coefficient Vq defined by: 



V = V -[ - ; - (2.16) 

V q J 1 — 2Vim*po 

This modified coefficient arises because of the change in the momentum density induced by the external 
field. For interactions with no momentum dependence, i.e., t\ = t2 = 0, we recover the results of 
Ref. @[|. 



Our formula for the response function ( 2.14| ) generalizes that of Garcia- Recio et al. [21] to which 
it reduces at zero temperature. 

B. Some properties of the Response Function 

A collective behavior will be observed when the response function exhibits a peak. Although this 
may be the case when the denominator in the response function goes through a minimum, the most 
familiar situation corresponds to the case where there is a pole in Eq. This occurs when the 

value loq of the frequency is such that 

1 = V U + 2VtU2 + tf(nl + IIoILO (2.17) 

The real part of the frequency cuq determines the energy of the resonance and its imaginary part the 
width. The corresponding relation is easily constructed when the imaginary part of uj is small so that 



a linearization of the previous equation near 9m(wo)=0 can be made 22]. In other cases a numerical 
construction of the strength function is necessary, as performed in some of the examples below. 

It is interesting to note that the symmetry energy coefficient a T of nuclear matter is, as expected, 
an important ingredient to determine whether or not there is a pole. Indeed one has the following 
relation between a T and the coefficients occuring in Eq. ( [2.17 ) 



Ur = 6^ + f \ V ° + TQ {t2(1 + X2) + 3tl(1 + + 2Fl ^) (2 ' 18) 

Note that this relation also involves the Fermi momentum kp = (37r 2 po/2) 1 ^ 3 and the momentum 
transfer q. 

1 A factor 2 is however missing in Eq. (13) of this reference, which implies that the interaction strengths in figs. 1 
and 2 of this reference must be divided by a factor 2. 
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The formula we have derived for the response function looks somewhat cumbersome and thus it 
seems worthwhile to explore whether approximations to the previous scheme can be derived. We have 
found that one such usefull scheme is provided by a Thomas- Fermi type appproximation in which one 
assumes that the change in the kinetic energy density r can be calculated from the proportionality 
relation r ~ p 5 ^ 3 , which leads to 

Sr(r, t)/r(r, t) = 55p(r, t)/3p(r, t). (2.19) 



Inserting this value into Eq. ( [2.11 ) one obtains the approximate expression for the response function 



„, \ Iln(ti;,q) 

U(u, q) = = 0V_1^ 2.20 

l-V U (oo,c L )-2V 1 n 2 (u J ,c l ) 



in which the coefficient Vq is defined by 



V = V - 2Vik 2 F . (2.21) 



This approximate formula was explored in Ref. It was found to reproduce correctly the main 

features of the exact formula derived in the previous section: existence of a zero sound at zero tem- 
perature with a disappearance of this mode at temperatures of a few MeV. 

In the limit u) —* and q —* 0, the response function becomes the static polarizability. In the 
isovector channel, the latter is related to the symmetry energy coefficient a T at temperature T. Indeed 
for q = uj = the energy density of nuclear matter in the presence of a small external field is given by 

H = H + ^^-(5p) 2 + e5p = H + —a 2 + ea. (2.22) 
Po Po 



By minimizing with respect to a we find 



a = , (2.23) 

yielding 

n(u, = 0,g = 0) = -^. (2.24) 

This formula can also be obtained by performing the limit u — > and q — > in the expression of the 
RPA response function ( 2.14| ), which provides a check of our formula for this quantity. 



The value of the strength distribution per unit volume S(uj) for the operator T3 exp(iq.r) is pro- 
portionnal to the imaginary part of the response function: 

%) = — 9mn(w,?). (2.25) 

7T 



It is also related to the photoabsorbtion strength distribution S a b s ||23| 



S abs (io) = — ^9mn(w, q). (2.26) 

7r i — e pw 



Furthermore, it satisfies the energy weighted sum rule |2 



oo 



-.2 



where 



du coS(u;) = -^p (1 + k) (2.27) 
2m* 



K = - 7 ^-p-(t 2 (l + 2x 2 )-t 1 (l + 2x 1 )) (2.28) 



is the enhancement factor arising from the the momentum dependent terms of the Skyrme interaction. 



III. Results and discussion 



The strength function is plotted for the Skyrme forces SGII f24|, SkM J2§, SIII > SI and SV H§ HI 
in Figs. 1, 2, 3, 4 and 5 respectively, as a function of the excitation energy to for various values of the 
temperature. The corresponding values of the Skyrme parameters are given in appendix A. In these 
figures we have chosen the value of the momentum transfer q in such a way that it corresponds to 



the dipole mode in lead-208 described by the Steinwedel- Jenssen model [28]. In this model neutrons 



oscillate against protons inside a sphere of radius R according to the formula 

(r\5p\r) = e sin(q.r) sin(atf), (3-1) 

where 

71" 

Taking R = 6.7 fm for lead-208 we find q = 0.23 fur 1 . 

We have checked in our numerical calculations that the energy weighted sum rule is well satisfied. 
Some examples are shown in Table 1, which compares the right hand side of the sum rule [c.f., Eq. 
( 2. 27] )] to the value of the integrated strength using Simpson's rule. It can be checked that there is an 



excellent agreement. There are some apparent exceptions however which correspond to the presence 
of a sharp zero sound peak which falls in between two meshpoints of the integration method. We have 
checked that the contribution of the pole just provides the missing strength in these cases. 

For interactions SGII, SkM and SIII the strength function shows a strong temperature dependence 
and exhibits a sharp peak at zero temperature. This peak occurs at a value of 18 MeV which is slightly 



Table 1: Right hand side (RHS) of the energy weighted sum rule (MeV x fm 3 ) compared to the 
integrated value of the strength m\ for T = 0, 3 and 6 MeV. 





RHS 


mi(T = 0) 


mi(T = 3) 


m\(T = 6) 


SGII 


52 


25 


52 


52 


SkM 


54 


15 


54 


54 


SIII 


49 


18 


49 


49 


SI 


43 


11 


43 


43 


sv 


121 


122 


121 


120 



higher than the observed one (we assume that we can rely on the Steinwedel- Jenssen model). One 
common feature of the three interactions SGII, SkM and SIII is that their effective masses are very 
close (respectively 0.78, 0.79 and 0.76, in units of the bare mass). 

For a Skyrme interaction with an effective mass closer to unity, such as SI (see Fig. 4), the 
resonance is found at a lower energy (15.5 MeV). For Skyrme forces with a very small effective mass, 
such as SV (m*/m = 0.36), the energy of the dipole resonance (see Fig. 5) is much higher (about 35 
MeV) than the observed value and the dipole strength is nearly independent of temperature. It still 
exhibits a weak collective behavior but for most of the energy range, the form of the response function 
is not very different from the imaginary part of the bare Lindhard function Ho, which is shown in 
Fig. 6. The imaginary part of IIo exhibits a nearly linear growth with energy which is also present 
in the imaginary part of the RPA response function, with moderate deviations occuring only in the 
resonance region. 

Let us now try to understand qualitatively the previous results about the position of the resonance 
and its evolution with temperature. The imaginary part of the bare response function IIo is known to 
show at zero temperatures two angular points at the following values of the energy (see appendix B 
and Fig. 6) 

qk F , q 2 ,„ Q s 

w ± = — ± 7T— ;■ (3.3) 

As temperature rises the corresponding discontinuities are smeared by the presence of the Fermi 
occupation numbers (Fig. 6). At these points the real part of the bare response function has a vertical 
slope (at zero temperature) and a maximum in between these points. For the particular value of the 
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momentum transfer we are considering (g=0.23 fm 1 ) and for a Fermi momentum kp =1.36 fm 1 one 
finds 

777 

w± = — x (15 ± 1.5)MeV. (3.4) 

m* 

Let us now assume for more simplicity (which is supported by numerical estimates) that the Lindhard 
function of order zero is the most important one for our discussion, so that the RPA response ( |2,14| ) 
can be approximated as 

n (w,q) 

II(a;,q)~ , 7. (3.5) 

V ,H) 1-Vbn (u;,q) 1 ; 

At this point we should note that the maximum value of the real part of the bare response is propor- 
tional to the effective mass 

max($fteIIo) = m*kpA(kp, q, T), (3-6) 

as can be seen from the equations given in appendix B. Since the maximum is smeared and reduced 
at high temperatures the function A decreases as temperature increases. In terms of this function the 
condition for the existence of a zero sound is 

1 < V m*k F A(T). (3.7) 

It turns out that for interactions with an effective mass close to unity and at zero temperature the 
condition (|3.7| ) is just satisfied. If this is the case, we expect the response function to have a maximum 
when the real part of the bare response is near its maximum, i.e., between the two points defined by 
Eq. ( |3.3|) . Since these points are quite close, a reasonable estimate of the resonance energy is 

rry 

0Jr^ — x 15MeV. (3.8) 

m* 

By comparing this formula with the results in Figs. 1-5, it can be seen to provide a good description of 
the resonance energies. Note that for other values of the momentum transfer q, a similar construction 
would give the following dispersion relation 

uj r ~ c x q, (3.9) 

with 

c ~ % = v F = 0.3 x c x ^L, (3.10) 

where we have taken kp=1.36 fm . This formula exhibits the important role played by the value of 
the effective mass. In actual calculations the sound velocity Co is slightly larger than the Fermi velocity 
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Table 2: Velocity Co of the collective mode compared to the Fermi velocity vf for the Skyrme forces 
SI, SGII and SV. 





SI 


SV 


SGII 


CO 


0.35 


0.77 


0.40 


V F 


0.30 


0.77 


0.35 



as can be seen from Table 2. The similarity between the two velocities shows that nuclear matter 
exhibits weaker collective effects than helium- 3. Indeed, in this last case the zero sound velocity is 
more than 3 times the Fermi velocity |j29|| . 

The dispersion relation ujr = coq implies that the resonance in light nuclei is located higher in 
energy according to the relation q = 2tt/R, i.e., evolves with mass number as ur ~ A -1 / 3 in agreement 
with the empirical formula |28|] . Since in our model the strength is concentrated in a single region this 
result also means weaker collectivity in light nuclei. 

Let us now show that the previous formulae also provide an explanation for the presence or not 
of a zero sound at zero temperature and for the different temperature dependences of the strength 
obtained for various interactions. Indeed from the relation between Vq, a T and m* we have 

m* V m* k 2 F 

ttPo = — a T ~ ( 3J1 ) 

ml m bm 

Taking a T =30 MeV and kp= 1-36 fm _1 this relation shows that the quantity m*Vo, which is the 
relevant one for our discussion, is much smaller for interactions with a small effective mass (0.4) than 
for interactions with m* /m ~ 1. Therefore if the condition fl3.7| ) is just satisfied at zero temperature 
for m*/m=l, such will not be the case for an interaction with m* /m=0A. Similarly, since the function 
A decreases with temperature it is also clear that the condition for the existence of a zero sound will 
eventually no longer hold at high enough temperatures (in actual calculations a few MeV). 

IV. Conclusion 

In conclusion we have found that the standard Skyrme forces SGII, SII and SkM give rise at zero 
temperature to a zero sound type collective mode exhibiting the usual dispersion relation E Tes = coq. 
The sound velocity for these forces is just slightly greater than the Fermi velocity, which implies a 
rather weak collective behaviour as compared to helium-3 for which there is a factor three between 
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the two velocities. For values of the momentum transfer corresponding in the Steinwedel-Jenssen 
model to the giant dipole mode in lead-208, we have found that the zero sound mode disappears at 
temperatures of a few MeV. This may be related to the saturation of photon multiplicities observed in 
some recent experiments ||. For Skyrme forces with a small value of the effective mass such as SV, we 
have found no zero sound at zero temperature and a weak variation of the strength with temperature. 
It is worthwhile noting that the previous forces all provide good descriptions of nuclear properties, 
such as binding energies and radii, all over the periodic table. In spite of this common property, they 
do give rather different predictions for the temperature evolution of giant resonances and also for 
their positions. Collective properties thus appear to give useful information on the effective nucleon- 
nucleon interaction. 

One limitation of our calculations is that although they do include the effect of the volume sym- 
metry energy, they ignore the effect of the surface symmetry energy wich is known to play a role, 
especially in light nuclei |3(J. They ignore as well shell effects which are also important in light nuclei. 
These effects indeed produce in this case a fragmentation of the strength in RPA calculations (see for 



instance the results of Sagawa and Bertsch [17] using the SGII force) whereas our model produces 
only a broad peak. Complete RPA calculations in finite nuclei (including heavy nuclei) would thus 
be of interest to complete the results of the present work. For such calculations we believe that the 
discussion we have presented would be a usefull guide. 
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Appendix A: Parameters of the interactions used 

We consider an effective Skyrme interaction of the following form: 



V12 = to(l + xoP*)S(n - r 2 ) + | (1 + Xl P a ) [6(n ~~ r 2 )k 2 + fc ,2 5(n - r 2 ) 

t 2 (1 + x 2 P„) k'.5(n - r 2 )k + |(1 + x 3 P (7 )p a 5(r 1 - r 2 ) 

o 



(A.l) 



where P a is the spin exchange operator. 

Denoting the single particle wave functions by (f>i(r,a,q), a and q being the labels for spin and 
isospin, the nucleon density Pg(r), the kinetic energy density r ? (r) and the momentum density j(r) 



11 



Table 3: Numerical values of the parameters to, h, t<i, t 3 , xq, x\, x 2 , X3 and a corresponding to the 
five Skyrme interactions SGII, SkM, SIII, SI and SV considered in this work. 





to 


h 


t 2 


*3 




Xi 


X 2 


S3 


a 




(MeV fm 3 ) 


( MeV fm 5 ) 


(MeV fm 5 ) 


(MeV fm 6 ) 












SGIP 


-2645. 


340. 


-41.9 


15595. 


.09 


-.588 


1.425 


.0604 


1/6 


SkM 6 


-2645. 


410. 


-135. 


15595. 


.09 











1/6 


SUP 


-1128.75 


395. 


-95.0 


14000. 


.45 








1 


1 


SP 


-1057.3 


235.9 


-100.0 


14463.5 


.56 








1 


1 


sv c 


-1248.29 


970.56 


107.22 


0. 


-.17 








1 


1 



a From Ref. j2§. b From Ref. [||], c From Refs. |26 



can be expressed as 

pg( r ) =Ei,a\Mr^,q)\ 2 

r q (r) =Ei,JV^(r,a,g)| 2 

j(r) = k (V - V) <P*(r',<J,q)Mr,^q)\r=r' 
The Hartree-Fock mean field Hamiltonian reads: 

where U q is a local potential given in the case of neutrons (q = n) by: 

U n (r,t) = t {(l + x i) p(r,t) - (xo + ±) p n (r,t)} 
+t 8 £{(2 + a)(l + f)p<^ 1 )(r l i))} 

-ts{(x 3 + i) (2p n (r, t)p a (r, t) + a(p 2 p (r, t) + p 2 n (r, t))p^ } 
+{^i(l + %) + 1*2(1 + f )} (r - i (V.j + j.V)) 
+{-|ti (si + i) + |t 2 (x 2 + |)} (r„ - £ (Vj, + j,.V)) 
+ {-i*i(2 + x i) + s(s + ^)} V 2 p + {|ti(i + x x ) + |i 2 (i + x 2 
with a similar expression for protons. The effective mass m* is given by: 



(A.2) 



(A.3) 



(A.4) 



(A.5) 
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Appendix B: Expression of the response function 

In this appendix we give the explicit expressions of the real and imaginary parts of the generalized 
Lindhard functions defined in Eqs. (|2.15|). For this purpose we need to define the following integrals: 



hN = _JL fey/?* (b.i) 

(2vr) 3 7 uj + zr] + e(k + q) - e(k) 

The imaginary parts of the generalized Lindhard functions are given by 

^n (^q)= -—^og l + e ^_ E+) 

Zmn 2 (u, q) = -™t (/V^log [lilt-El] + ^ C 1 + ^~ E+) ) ~ ^ C 1 + e^)) 



Qrrtll^a;, q) = — qyj2m*E + ^smI\.2 + 25m/4 — 2y/2m* E + qQml2 



where Li2 is the Euler dilogarithmic function [31] 



and 



2g 2 v 2m* 

The expressions of the functions Q171I2 and Qml^ are : 

m* 3 / 1 -I- pP(^~ E- ) 

Zml 2 (u, q) = -— 2 [PE + log i ± ^_ E+) + Li 2 (l + e^ + )) - Li 2 (l + e^-)) 

2m* 4 9 ( l -|- e £0--E-) roc 1 _|_ e -p(zE+-u-n) ^ 

9m/ 4 (u;, g) = log — — , . + 2 / dz z log ■ 



?rg/? + ^ ta 1 + eP0*-E+) Ji 1 + e -^ zE +~^ 

At zero temperature the real parts of the H 2 N are given by 

M (T = 0)= ^(-l + ^[^ + ) + 0(x_)]) 
M 2 (T = 0)= !^Ml(-3 + x + x_ + xi + x?.+ 



2q 



q 2 m*E + ^en (T = 0) + ^m*q 2 k 3 F 



(B.2) 



LioCO / ^yd/ (B.3) 



(B.5) 



[(1 - x\ - 2x + x„)(f>(x + ) + (1 - xt - 2x + x_ )</>(>-)] j (B.6) 

M 4 (T = 0) = 2 (Ke/ 4 (T = 0) - 2q^2m*E + $teI 2 (T = 0) + 

1 
3^2 
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(1 + x\)^{x + ) + (1 + x± + 



(B.7) 



where 

Ke / 2 (T = 0)= ^^-S-x+x.-x^ + x 2 

$Re/ 4 (T = 0) = ( (1(1 + xl + xt) + ^(1 + x 2 ) + 

1, 9 4 . 5g 1,1 * 1^ 2m*wi_ 8g 2 x + , 

-(1 + xi + xiU(x+) + — ^- + -(-x^ + ~xl H p5 h , )+ 

6 V + +; v +; 3/c F 3 V 3 _ 3 + fcf, /c 2 ; 

x^ x+ 2m*u;xi 4m* 2 u; 2 x_\ 
3 3 J 

, . q m*uj , . 2 x — 1 

with x± = — — ± — - — , and <p(x) = (1 — x ) log . 

2k F qk F x + 1 

For non zero temperature the expression of 5ftell 2 jv (N=0,l,2) is an average of the zero temperature 
functions calculated for the same values of u> and q, but with various values of the Fermi momentum 
k F distributed with a weight factor which is just the derivative of the Fermi occupation number. 
Explicitly one has the following formula: 

Ken (w, q,T) = -J Ken (w, q,T = 0,k F = k) df(k, T) 

5?en 2 (w, q, T) = - J M 2 (w, q,T = 0,k F = k) df(k, T) (B.8) 
M 4 (w, q,T)=-J M 4 (w, q,T = 0,k F = k) df(k, T), 

where f(k,T) is the occupation number. For the case of zero temperature we have 

df(k) = - 5(k-k F ) dk, 

yielding the above expressions for these functions. 

We would like also to show the relation between our definition of generalized Lindhard functions 
( II 2 j ) and those of Garcia- Recio et al. plj (n 2 jv) defined in the limit T = 



n = 4n 

n 2 = 4n 2 - 2g 2 n (b.9) 

n 4 = 4fi 4 - 4g 2 fi 2 - q 2 m*p + q 4 U Q + (2m*ujq) 2 Uo 



Our definition of Vq and V\ and their W\ and W2 are related by 

A 8 (B.10) 
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Figure captions 

Figure 1 Distribution of strength per unit volume for the operator exp(zq.r) (in fm~ 2 ) as a function 

of the energy to (in MeV) for a momentum q = 0.23 fm _1 and for different values of the temperature 

T = — > 6 MeV, in the case of the interaction SGII. 

Figure 2 Same as Fig. 1 for the force SkM. 

Figure 3 Same as Fig. 1 for the force SIII. 

Figure 4 Same as Fig. 1 for the force SI. 

Figure 5 Same as Fig. 1 for the force SV. 

Figure 6 Imaginary part of Ho ( in fm~ 2 ) as a function of the energy uj ( in MeV), in the case of SV 
force, for q = 0.23 fm -1 ; the temperatures are from T=0 to T=6 MeV. 
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